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Motivated by the partial differential equations of mixed type that arise in the 
reduction of the Einstein equations by a helical Killing vector field, we consider 
a boundary value problem for the helically-reduced wave equation with an arbi- 
trary source in 2+1 dimensional Minkowski spacetime. The reduced equation is 
a second-order partial differential equation which is elliptic inside a disk and hy- 
perbolic outside the disk. We show that the reduced equation can be cast into 
symmetric-positive form. Using results from the theory of symmetric-positive dif- 
ferential equations, we show that this form of the helically-reduced wave equation 
admits unique, strong solutions for a class of boundary conditions which include 
Sommerfeld conditions at the outer boundary. 



1. Introduction 

Physical systems are typically governed by partial differential equations (PDEs) of a 
fixed type: elliptic, hyperbolic, or parabolic. The mathematical properties of such equa- 
tions have been extensively investigated (see, e.g., Refs. [1,2]). Considerably less is known 
about PDEs of mixed type, by which we mean equations whose type is different in different 
subdomains of the domain of interest, e.g., elliptic in one region and hyperbolic in another 
[3]. Compared to elliptic, hyperbolic or parabolic equations, mixed type equations are 
rather unusual, both in the boundary conditions that can be imposed to get existence and 
uniqueness of solutions as well as in the regularity of solutions that are obtained. More- 
over, the lower-order terms in equations of mixed type take on a more significant role than 
in equations of fixed type. This latter feature means that it is difficult to obtain general 
results about PDEs of mixed type; to a large extent, one must investigate each set of 
equations, each set of boundary conditions, etc. separately. 

In relativistic field theory on a fixed spacetime, mixed type equations occur after 
performing a symmetry reduction of hyperbolic PDEs with respect to an isometry group 
which has an infinitesimal generator that changes type from timelike to spacelike. In 
generally covariant theories such symmetry reductions may yield PDEs of mixed type in 
appropriate gauges. An important example of the latter type, currently of considerable 
interest in gravitational physics, arises in the quasi-stationary approximation to the 2- 
body problem in general relativity [4,5,6]. There one is interested in solving the Einstein 
equations for spacetimes admitting a helical Killing vector field. The helical Killing vector 
field, which represents a rotating reference frame, will be timelike near the bodies, and 
spacelike far from the bodies. The reduced Einstein equations (modulo gauge) can be 
expected to include non-linear PDEs of mixed type on the 2+1 dimensional manifold 
of orbits of the Killing vector field. The reduced equations must be solved numerically, 
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but one naturally desires as much a priori information about existence and uniqueness 
of solutions, regularity of solutions, admissible boundary conditions, etc. as one can get. 
The quasi-stationary approximation to the relativistic 2-body problem is the principal 
motivation behind the work presented here. 

As a warm-up for numerically solving the Einstein equations in the quasi-stationary 
approximation, the wave equation and some of its non-linear extensions, reduced by the 
assumption of helical symmetry, have been examined both analytically and numerically in 
[7] . The helically- reduced wave operator is elliptic inside a "light cylinder" and hyperbolic 
outside the cylinder. It was shown in [7] that one can give a formal series solution of 
the helically-reduced wave equation for a source consisting of a pair of equal and opposite 
point charges, placed symmetrically with respect to the axis of helical symmetry. These 
analytical solutions satisfy Sommerfeld conditions at an outer boundary (which may be at 
infinity) [8] . The choice of such boundary conditions can be motivated on physical grounds, 
and the apparent analytic existence of unique solutions in the linear case and numerical 
solutions in the non-linear case gives confidence that the helically-reduced equation can 
be treated as a boundary value problem. However, from a mathematical point of view it 
is not immediately clear a priori why such boundary conditions are admissible, that is, 
why one should expect unique solutions to exist. Our goal here is to understand existence, 
uniqueness, regularity, etc. — in short, the well-posed nature of this problem — from a 
general point of view that does not rely upon explicitly constructing a solution to the PDE. 
The idea is that such a point of view can be used to better understand the helically-reduced 
Einstein equations, which will not yield so easily to a direct assault. 

In the past, certain boundary value problems of mixed type have been addressed us- 
ing the theory of "symmetric-positive" differential equations [9], which can be viewed as 
a generalization of elliptic and hyperbolic equations. Friedrichs [9] and Lax and Phillips 
[10] have given the basic existence and uniqueness results for linear symmetric-positive 
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equations. In this note we shall show that the helically-reduced wave equation with arbi- 
trary sources in 2+1 dimensions can be cast into symmetric-positive form. We can then 
deduce existence and uniqueness results for a class of boundary conditions that include 
the Sommerfeld conditions used in [7]. These results provide support for the proposition 
that boundary value problems of mixed type arising from helical symmetry reductions - 
such as arise in the relativistic two-body problem — are well-posed. 

2. The helically-reduced wave equation 

We consider 2+1 dimensional Minkowski spacetime (R 3 , rf) and a helical Killing vector 
field £. There will exist an inertial-Cartesian coordinate chart (t, x, y), such that 

rj = —dt <E> dt + dx <8> dx + dy <E> dy, (2.1) 

and 

for some constant O. In the corresponding inertial-polar coordinates (t, r, <fi) we have 

77 = -dt ®dt + dr®dr + r 2 d<j) <g> #, (2.3) 

and 



Note that ^ is not of a fixed type: 



> for r > Yi 



V (C,0=r 2 n 2 -l {=0 forr=l (2.5) 



n 

< for r < yj. 



We will call the set of spacetime events with r = ^ the cylinder. 
The wave equation with source j: R 3 — > R is given by 



n* = j. (2.6) 
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In the inertial-polar chart, the wave operator acting on a function <£>: R 3 — > R takes the 
form 

□$ = -«9 2 $ + ^dr(rdr^) + (2.7) 

We now restrict attention to fields and sources which are invariant under the 1- 
parameter isometry group G generated by £. This is equivalent to requiring 

m = o=m- (2.8) 

These conditions imply that $ and j define functions on the manifold of orbits R 3 /G, which 
shall be denoted by \1/ and /, respectively. Because the source j is assumed G-invariant, 
and because G is an isometry group for the spacetime, the wave equation (2.6) admits G 
as a symmetry group and defines a PDE relating \1> to / on R 3 /G [11,12]. To obtain this 
differential equation in local coordinates we proceed as follows. Group invariants on R 3 
are functions of r and 

(p:=(j)-Qt, (2.9) 
which define polar coordinates on R 3 /G ~ R 2 . In particular, granted (2.8), we have 

$(t, r, 4>) = *(r, <p), j(t, r, <j>) = f(r, ip). (2.10) 

The reduced field equation on R 3 /G can be obtained by substituting (2.10) into (2.6), 
which gives 

ldr(rd r y) + ^x(r)dfo = f, (2.11) 

where 

X(r) = l-0 2 r 2 . (2.12) 

Note that the light cylinder on R 3 projects to a light circle at r = ^ on R 3 /G. 
Evidently, (2.11) is elliptic inside the light circle and hyperbolic outside the light circle, 
which is a consequence of the changing character (2.5) of the Killing vector £. Thus (2.11) 
is a PDE of mixed type. 
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In [7] (2.11) is solved on a disk of radius R with source / corresponding to two equal and 
opposite "scalar point charges" placed symmetrically relative to the origin. Sommerfeld 
conditions are imposed at the boundary of the disk; and it is required that ^ vanishes at 
the origin. The solution is given as a formal infinite series. As noted in [7], despite the 
mixed type of the PDE and, in particular, despite the fact that a Sommerfeld condition 
was enforced at an outer boundary in the region where the PDE is hyperbolic, a unique 
solution exists. Here we provide a somewhat more general version of this result using the 
theory of symmetric positive equations. We consider (2.11) on a domain [13] e < r < R, 
< e < ^. We allow for a general source / and we employ a class of boundary conditions 
that include Sommerfeld conditions at r = R, such as considered in [7]. Specifically, we 
will impose the following boundary conditions: 

*(e,^) = 0, T((p)Rd r ^(R,(p) + <r((p)d^(R,(p) = 0, ar ± (2.13) 

where a and r represent smooth functions on the outer boundary r = R. Sommerfeld 
boundary conditions correspond to setting t = 1/R and a = ±0. 

Although we explicitly consider homogeneous boundary conditions (2.13), because we 
allow for an arbitrary source / in (2.11) a large class of inhomogeneous boundary conditions 

#(e, ip) = k(<p), r(ip) Rd r ^(R, ip) + a(<p) d^(R, ip) = l(<p), (2.14) 

can also be accommodated. This is done by choosing a smooth function A = A(r, ip) which 
satisfies the inhomogeneous boundary conditions (2.14) and then redefining 

*^* = *-A, f^f = f- l -d r (rd r K) + ^x(r)d*A. (2.15) 

^ now satisfies (2.11) with source / and homogeneous boundary conditions (2.13), to which 
our results apply. 

3. Symmetric-positive PDEs 

Existence and uniqueness results can be obtained for linear PDEs of mixed type if they 
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can be cast into first-order, symmetric-positive form with appropriate boundary conditions 
[9,10]. Here we summarize the results from [9,10] which we shall need. 

For our purposes, the data used to define a symmetric positive system of equations 
with admissible boundary conditions will be taken to be [14]: 

(i) A smooth manifold M with smooth boundary <9M; we set M = M U dM. 

(ii) A smooth scalar density of weight-1 on M, denoted by uj. 

(iii) A finite-dimensional real vector space V with scalar product (•,•)• 

Let u: M — > V. We consider a first-order system of differential equations for u on M 
of the form 

Lu = A a V a u + Bu = h, (3.1) 
where h: M — > V and, at each x G M, A a and B are linear transformations: 

A a (x):T*M xV ^V, B(x):V^V. (3.2) 

For simplicity we assume that A a and B depend smoothly on x E M. The differential 
operator V a is the exterior derivative on functions defined by any basis for V. 

Definition 3.1. The system (3.1) is symmetric-positive if: (1) A a defines a symmetric 
operator (with respect to the scalar product on V): 

(A a (x)v a )* = A a (x)v a , G T*M and V x G M, (3.3) 

and (2) the linear operator K(x): V — > V, defined by 

K = B-^V a (uA a ), (3.4) 

has a positive-definite symmetric part: 

K(x) + K*(x) > 0, \/x G M. (3.5) 
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We remark that V a in (3.4) is defined by its unique torsion-free extension to vector densities 
of weight-1 on M taking values in V <E> V*. 

A class of boundary conditions on u:M ^ V have been determined such that there 
exist unique solutions to symmetric-positive systems of PDEs. Following Friedrichs, we 
call these boundary conditions admissible. They are defined as follows. 

Definition 3.2. Boundary conditions on u: M — > V are the requirement that, at each 
x G dM, u(x) takes values in a linear subspace N(x) C V, which varies smoothly with x. 

We fix an outwardly oriented normal 1-form, n a , to the boundary dM. This 1-form is 
uniquely determined up to multiplication by a smooth positive function on the boundary. 
We define [15] 



P = n a A a 



(3.6) 



Definition 3.3. Let u: M — > V be subject to the system of equations (3.1). The boundary 
conditions u(x) G N(x) on dM are admissible if N(x) is a maximal subspace such that 
the quadratic form u — > (u, (3u) is non-negative at each x G dM. 

Note that the admissibility of a set of boundary conditions does not depend upon the 
specific choice of outward covariant normal to the boundary. In the sequel we will make 
use of the following convenient characterization of admissible boundary conditions, which 
is due to Friedrichs [9]. 

Proposition 3.4. Admissible boundary conditions, u(x) G N(x) on dM are equivalent 
to the linear boundary conditions {5<i u — on dM, where (32 arises from a decomposition 

P = Pl+p2 (3-7) 
such, that for all x G dM, (i) every v G V can be decomposed via 

v = v\+ V2, P\V2 = foyi = 0, (3.8) 
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and (ii) 



(3.9) 



has a non-negative symmetric part: 



/I + > 0. 



(3.10) 



We now summarize the existence and uniqueness results of [9,10] for symmetric-positive 
systems. We say that a mapping u: M — > V is in L 2 (M, V) if 



Definition 3.5. Let the mappings u: M -> V andh:M be in L 2 (M, V) ; u is a strong 
solution to (3.1) satisfying the boundary conditions u(x) G N(x) on dM if there exists a 
sequence of functions G C°°(M, V), satisfying the boundary conditions u^{x) G N(x) 
on dM, such that 

Ufc — > u, and Z/u& — > /i 

in the L 2 (semi-)norm (3.11). 

Theorem 3.6 (Friedrichs, Lax & Phillips). If (3.1) is symmetric-positive then it 
admits a unique, strong solution satisfying admissible boundary conditions. 

We remark that the theorems appearing in [9,10] use (piecewise) continuously differen- 
tial) le functions to define strong solutions. However it is straightforward to check that 
the relevant results go through for {u^} G C°°(M, V), which we use here. 

4. A symmetric-positive system for the helically-reduced wave equation 




(3.11) 



Here we show that the helically-reduced wave equation can be expressed in symmetric- 
positive form. We choose M to be an annulus, M = {(r, ip)\e < r < R}, equipped with the 
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metric 

g = dr ® dr + r dip dip, (4-1) 

and associated density: 

uj = (det g)^ = r. (4.2) 
We set V = R and equip it with the standard scalar product 

(it, v) = u\v\ + U2V2- (4.3) 

We consider the following first-order system 

1 1 

-d r u 2 + ~2X<Vn = /, (4-4) 
1 1 

-d r U\ - ^2#yjU2 = 0. (4.5) 

Setting 

ui = d^, ix 2 = rd r V, (4.6) 

all classical (C 2 ) solutions to (2.11), (2.13) are solutions to (4.4), (4.5). The solution it 
thus obtained satisfies the boundary conditions 

Ul (e,ip)=0, a((p) Ul (R,(p) + r(ip)u 2 (R,ip) = 0, ar ^ 0. (4.7) 

Conversely, given a classical (C 1 ) solution to (4.4), (4.5) satisfying boundary conditions 
(4.7), the function * defined by 

f r 1 

tf(r >¥ >) = J dr'-jU 2 {r',cp), (4.8) 

satisfies (4.6) and hence (2.11), (2.13). In this sense the equations (2.11), (2.13) are 
equivalent to (4.4), (4.5), (4.7). We write the system (4.4), (4.5) as 

(A a V a u -~h)=0, (4.9) 
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where 



and 



«=(SJ- (411) 

Now consider the following first-order system, 

L(A a V a u - h) = 0, (4.12) 

where 

L=(l -«), (4.13) 
and a = a(r) and c = c(r) are smooth functions on M. We require that a and c satisfy 

a 2 + c 2 x^0. (4.14) 

so that (4.12) is equivalent to (4.9). We write (4.12) as 

A a V a u = h, (4.15) 



where 



and 



A r = -(-" x ;Y ^=4r* ?J, (4.16) 



r \ a c J ' r 2 V C X ~~ a 



We now show that the functions a and c can be chosen so that the system (4.15) is 
symmetric-positive. The linear transformations defined by A a are symmetric for any choice 
of a and c. Using (3.4) we find that 



- JV* - I ( dr(cx) -0, 
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Necessary and sufficient conditions for (3.5) are 



d r C < 0, (drC)dr(cx) < -(d r af. 



(4.19) 



A pair of smooth functions that satisfy (4.19) and (4.14) for < e < r < R are of the form 



provided the constants \a\ and a > are chosen large enough, as is easily verified from 



Thus, with these choices for a and c, the equations (4.15) are symmetric-positive. We 
summarize the preceding discussion as follows. 

Proposition 4.1. The equations (2.11), (2.13) are equivalent to the symmetric-positive 
system (4.15), (4.7), (4.20). 

5. Admissible boundary conditions 

We now consider the boundary conditions (4.7) for which we have the following result. 

Proposition 5.1. The constants a and a in (4.13) and (4.20) can be chosen so that the 
boundary conditions (4.7) for (4.15) are admissible. 



a = const, c = —a + e 



(4.20) 





(4.21) 



Proof: 



For the normal 1-form to the boundary we use 




Rdr at r = R 
—edr at r = e. 



(5.1) 



We then have 




(5.2) 
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We consider boundary conditions of the form 



au\ + TU2 = on dM, (5.3) 

where 

cr = l, r = 0, at r = e, (5.4) 

and 

a = a((p) ^ 0, t = r((p) ^ 0, at r = R. (5.5) 
Following Proposition 3.4, we set 

P 1 = ±N ( - T l CX - ara ^c X + * 2 a\ (52 = ±N (-o 2 c X + ara -arc X +fa\ 

\ t a — arc —ara + a c ) \ a a + arc ara + r e J 

(5.6) 

where 
and 

= 01+02- (5-8) 

In (5.6) the plus/minus sign is to be used at the outer/inner boundary. The boundary 
conditions (5.3) are equivalent to 

2 u = 0, on dM. (5.9) 

At the inner boundary we have 

i(M + ,*) = (- c(£ « £) _ c ° (e) ). (5.10) 

We choose a sufficiently large so that c(e) < and we assume the inner boundary is within 
the light circle so we have %(e) > 0. Therefore {ji + /j,*) is non-negative at the inner 
boundary; the inner boundary condition is admissible. At the outer boundary we have 

-(^ + ^)=ivf- 2(Tra+((T 2 2 2 r2)CX -"Tiff* 2, ) ■ (5-H) 

2^ P ' V -arOVc -2ara + (a 2 - r 2 )c J r=R V ; 
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Evidently, the outer boundary conditions are admissible in the case or > provided a is 
chosen sufficiently negative [16]. If ar < the boundary conditions are admissible provided 
a is chosen sufficiently positive [16]. Thus the outer boundary conditions are admissible 
provided or ^ there. 

We remark that the admissibility of the outer boundary conditions did not depend 
upon the location of the outer boundary. Therefore the outer boundary conditions can 
be imposed outside, inside, and even on the light circle. We also note that the proof of 
Proposition 5.1 shows that at the outer boundary neither Dirichlet conditions (a = 1, r = 
0) nor Neumann conditions (a = 0, r = 1) for \1/ are admissible. 

6. Conclusions 

From Propositions 4.1, 5.1 and Theorem 3.6 we have the following basic existence and 
uniqueness result for the helically-reduced wave equation in its symmetric positive form. 

Theorem 6.1. The system (4.15) admits a unique, strong solution satisfying the boundary 
conditions (4.7). 

While this theorem only establishes existence of distributional solutions to the first- 
order form of (2.11), (2.13), it does imply that classical solutions to (2.11), (2.13) are 
unique. Considerable additional analysis is needed to establish existence of C 2 solutions 
to (2.11), (2.13). However, using Theorem 6.1 the following regularity properties of \1/ can 
be immediately inferred. 

Let us define H l (M) as the completion in the if 1 norm, 




(6-1) 



of the space of smooth functions \1/ satisfying the boundary conditions (2.13). 



We can then deduce the following from Theorem 6.1. 
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Corollary 6.2. Let u be the strong solution to (4.15), (4.7). The function ^, defined by 

* = J* dr' y,u 2 {r') = I(u), (6.2) 
is in H l (M) with (distributional) derivatives given by 

ui = d^, u 2 = rd r ^. (6.3) 

Proof: 

I is easily verified to be a bounded linear transformation from the dense subspace 
C°°(M, V) C L 2 (M, V) to L 2 (M), so / can be extended to all of L 2 (M, V) thus defining 
\1/ G L 2 (M) via (6.2). Let u n : M — > V be the sequence of smooth maps that converges to 
the strong solution u. Clearly, 

I(u 2n ) = 0, (6.4) 
r=e 

so that the inner boundary condition is satisfied. We have 

drl{u 2n ) = -u 2n -> -«2 e L 2 (M, V), (6.5) 
r r 

and it is straightforward to verify that 

d^I(u 2 n) = I(d^u 2n ) ^«iG L 2 (M, V). (6.6) 

Thus \1/ and its first derivatives are in L 2 and satisfy 

ui = <9^, u 2 = rd r ^, (6.7) 

so that the outer boundary conditions are satisfied as well. 
□ 



Physically, the source / in (2.11) cannot be known with perfect precision. Furthermore, 
one may only have an approximately correct source appearing in a numerical solution. It 
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is therefore important to note that the solution to (4.15), (4.7) depends continuously upon 
the source h, so small changes/errors in the choice of h lead to correspondingly small 
changes in the solution u (or \1/). (Here "small" is defined by the L? norm.) To see this, 
we define a linear mapping S: L 2 (M, V) — > L 2 (M, V) that associates a solution u = S(h) 
of (4.15), (4.7) to each source h. The existence of this mapping follows from Theorem 6.1. 
The uniqueness of strong solutions to (4.15), (4.7) implies that the mapping S is closed 
and hence continuous by the closed graph theorem (see e.g., [17]). From Corollary 6.2 the 
mapping I:L?{M, V) — > H^(M) is bounded - hence continuous - and we then have the 
following corollary. 

Corollary 4.3. ^ G fl" 1 (M) defined by (4.8) depends continuously upon the source f G 
L 2 (M). 
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